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Abstract. Wc study the Einstein relation between diffusion and response to an 
external field in systems showing supcrdiffusion. In particular, wc investigate a 
continuous time Levy walk where the velocity remains constant for a time r with 
distribution Pt{t) ~ . At varying g the diffusion can be standard or anomalous; 
in spite of this, if in the unperturbed system a current is absent, the Einstein relation 
holds. In the case where a current is present the scenario is more complicated and the 
usual Einstein relation fails. This suggests that the main ingredient for the breaking 
of the Einstein relation is not the anomalous diffusion but the presence of a mean drift 
(current). 
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1. Introduction 

The fluctuation-dissipation relation (FDR) is one of the fundamental results of the 
statistical physics [1]. In his celebrated work on the Brownian motion, Einstein gave 
the first example of FDR. In the absence of external forcing, at large times, one has 

= and {x^{t)) ~ 2Dt, (1) 

where x is the position of the colloidal particle, D is the diffusion coefficient and the 
average refers to the unperturbed dynamics. If a small constant external force S is 
applied, one obtains a linear drift 

{x{t))e ^ fi£t, (2) 

where {■ ■ ■)e is the average on the perturbed system, and ^ is the mobility. The 
remarkable result obtained by Einstein is the proportionality between {x(t))s and the 
mean square displacement (MSD) {x'^{t)): 

{x{t))e ^ ^ 

namely /i = /3-D, where /3 is the inverse temperature and we have taken the Boltzmann 
constant /c^ = 1. 

In the last decades many researches have been devoted to the study of anomalous 
diffusion, where, instead of ([T]), one has 

{x\t)) ~ t^'' with u ^ 1/2, (4) 

see for instance [21 [3], HJ |5l [6l [7] . The case z/ < 1/2 is called subdiffusion while if z/ > 1/2 
we are in the presence of super diffusion. It is quite natural to wonder about the validity 
of the FDR in these anomalous situations. Important results showing the proportionality 
between {x{t))£ and {x'^{t)), which we refer to as the Einstein relation "at equilibrium", 
have already been obtained for both superdiffusive [SI El C] and subdiffusive [21 (TUl E] 
anomalous dynamics. Differently, the situation where the drift is compared to the mean 
square displacement in a state which is already "out of equilibrium", either due to a 
current [11] or due to dissipation [T2], has been studied only in the subdiffusive cases. 

The aim of this letter is to discuss the validity of the Einstein relation in equilibrium 
and out-of-equilibrium situations in the presence of superdiffusive dynamics. In 
particular, we consider a Levy walk collision process [13] and we show that the Einstein 
relation is violated when the perturbation is applied on a reference state where a current 
is already present. 

2. The model 

We consider an ensemble of probe particles of mass m endowed with scalar velocity v 
and position x. Each probe particle only interacts with particles of mass M and velocity 
V extracted from an equilibrium bath at temperature T. We assume that the scattering 
probability does not depend on the relative velocity between the probe particle and the 
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colliders, as, for instance, in the case of Maxwell-molecule models [2]. Velocity of the 
probe particle changes from v to v' at each collision, according to the rule: 



where 7 = (C — a)/(l + C); with ( = m/M, and a is the coefficient of restitution (a = 1 
for an elastic collision). The velocity V of the bath particles is a random variable 
generated from a Gaussian distribution with zero mean and variance T/M: 



The elementary step of the dynamics is made by: i) a flight, x(t+r) = x(t)+w'r+l/2£^r^, 
where x{t) is the position of the probe particle at time t, with r taken from a distribution 
Pt{t) and £ a constant acceleration, followed by ii) a collision v' = '-fv + {1 — '-f)V , with 
V taken from the Gaussian distribution (jS]). In the speciflc case a = 1 and M = m, one 
has 7 = and the collision rule ([5]) results in a random update of the velocity according 
to the distribution ([6]). The duration of each flight, r, is an independent identically 
distributed random variable with probability 



with g > 1. This kind of process is called Levy walk collision process [13], and may 
be interpreted as due to scattering centers randomly distributed on a fractal spatial 
structure, as for instance in the case of molecular diffusion in porous media |T5]. If 
a 7^ 1 or m 7^ M, with a ^ (, a. dependence on the last velocity before the collision 
remains. In this case velocity correlations can be measured in the system, as discussed 
in Section [3711 

According to the dynamic rules of the process described above the displacement 
of the probe particle is always flnite in a flnite time. The anomalous dynamics of such 
a model has been studied in [16], showing that the process is an example of "strong" 
anomalous diffusion, namely that it is not possible to flnd a scaling for the probability 
density function (PDF) of displacements. Such a collision process becomes a standard 
diffusive system when P^- decays fast enough: in this regime the dynamics is qualitatively 
equivalent to that of a system with exponential Pr, studied for instance in [171 [HI [19] . 
Levy walk collision processes have been thoroughly studied, see for instance [7], where 
the behavior of the mean square displacement has been obtained analytically. 

We recall here a simple argument already presented in [16] to study the asymptotic 
behavior of higher order moments of displacement, and that can be easily applied also to 
the case with an external perturbing fleld that we are going to discuss here in Section [3] 
In order to obtain in a simple way the dominant asymptotic behavior of {x^{t)), we 
introduce a cut-off t^ for Pr{T~): 



= 7f + (1 — 7)V, 



(5) 




(6) 




-g 



(7) 




(8) 
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Assuming x = as initial condition for each trajectory, the mean square displacement 
after the time t, where N{t) collisions occurred, can be written in full generality as: 



(At)) 



N{t) 

ViTi 

1=1 



2- N{t) N{t) 



^EK^'>+2^WE(w.ro). (9) 

1=1 1=1 

Here, Vi denotes the velocity of the probe particle after the i-th collision, Tj is the time 
elapsed between the collisions i and i + 1 and N{t) is the average number of collisions 
occurred up to time t. The average (■ ■ ■) is taken over the distributions ([6]) and ([7]). 
From Eq. ([8]) we have, for n + 1 — (7 > 0, 

{r^c ~ tr'^^ (10) 
where (■ ■ ■)c denotes an average over the distribution ([8]) with the cut-off tc- 

We start by considering the case with independent velocities Vi, corresponding to 
a choice of parameters such that 7 = 0. Then, estimating (x'^it)) at a time t ^ t^, 
so that the average number of collisions along the trajectory can be approximated to 
N{t) ^ t/ {r)c, and considering that the cross terms in Eq. ([9]) are zero, we can write 

{x\t))^^{v'){r')^. (11) 

In the case g > 3, both (r)c and (t^)c are finite, even in the limit tc — ^ 00, so that we 
find the simple diffusive behavior {x^{t)) ~ t. For t <tc and 1 < g < 3, instead of (fTTjl 
we expect 

{At))-t'r (12) 

One can easily find the exponent u with a matching argument. Comparing f|TT]) and f|T2l) 
at t ~ tc and using (11 01) . we obtain i/ = 1/2 for (7 > 3, z/ = 2 — g/2 for 2 < (7 < 3, and 
z/ = lforl<(7<2 (logarithmic corrections appear at the values g = 2 and g = 3 [7]). 



3. Einstein relation 



The main concern of our study is the question about the validity of the Einstein relation 
for superdiffusive anomalous dynamics. In particular such a relation can be checked in 
two different situations: 

(A) The drift due to the external force is compared with the MSD of the probe 
particle in the absence of any pulling force. This case corresponds to a fluctuation- 
dissipation experiment realized by switching on from zero the external perturbation. 

(B) The drift is compared to the MSD in a state where a current is already present. This 
procedure corresponds to increase the intensity of the perturbation in a state already 
perturbed and compare the average current in this state with the fluctuations in 
the initial reference state. 

In the following, we will refer to situation (A) as a test of the fluctuation-dissipation 
relation at equilibrium while to case (B) as a test out of equilibrium. We will show that 
these two cases are very different. 
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Figure 1. Log-Log plot of the MSD in the case of superdiffusive dynamics with 
g = 3/2, {x'^{t)) (black line), and of the drift {x{t))£ (red dashed line) in the case 
of uniform external field £ = 1.5. Both quantities behave as ~ and the Einstein 
relation is verified. 



3.1. Perturbation of a state without current 

The argument used in Sec. |2]to study the MSD can be apphed to the drift, yielding 



Nit) 



i=l 



t 



{r)c 



£ t 



{r')c. (13) 



which perfectly matches the result for the MSD found in Eq. (1 11 1) . Therefore, when we 
perturb an equilibrium state, namely a state without currents, we have for any value of 

TTTrf = const. 14 

Let us note that the Einstein relation holds quite generally, namely it persists even if 
we make our process non-Markovian by allowing some memory across collisions, that is, 
we put 7 > in the collision rule Eq. ([S]), preventing a complete reshuffling of velocities. 
In this case the MSD reads as: 



Nit) 



t 



h3 



Nit) 



(15) 



because the collision times are still not dependent on the velocity but correlations 
between the velocities must be taken into account. Note that, even in the presence 
of non-independent {vi}, if '^^{viVo) yields a finite contribution at large times the 
second term on the right of Eq. (fT5|) grows as (t)^, namely it is subdominant compared 
to (t^)c in all situations. If the exponent of the flight time distribution is 1 < < 2, 
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Figure 2. Log-Log plot of the MSD (black line) in the presence of a constant external 
field E = 1.5 and drift (red dashed line) with the same value of the field, in the case of 
superdiffusive dynamics with g = 5/2. We observe a breaking of the Einstein relation 
with leading asymptotic behaviors, l\6x{tyf)^ ^ t^/"^, see Eq. p^. and {x{t))£ ^ t^/'^, 
see Eq. 

then (r^)c ~ t^~^ > ~ (r)^; if 2 < (7 < 3 we have that (r)^ is a finite number, 

while for g > 3 we recover a simple diffusive behavior for both the MSD and the drift. 
In particular, in the case g > 3 the presence of correlations of velocities amounts to a 
renormalization of the diffusion coefficient. 

We conclude that the asymptotic behavior of the MSD is the same as in Eq. ( fTTl) . 
From Eq. (fT3|) one can see that no cross correlations between velocities at different times 
appear, so that the drift is not influenced by velocities correlations across collisions. 
Therefore, we can conclude that the Einstein relation is preserved for all values of the 
exponent g of the flight times distribution. This example shows how the Einstein relation 
is a stable result, valid under quite general conditions also in the presence of anomalous 
dynamics [71[20l[II]. 

3.2. Perturbation of a state with a current 

In order to study the Einstein relation out of equilibrium, let us flrst consider a simple 
Gaussian process, namely the Brownian motion of a colloidal particle when we add 
a constant force pulling the particle. In this case it is sufficient to replace the MSD 
around the average position {[Sxit)]'^)^ = {x^it))^ — (x(t))| for the simple MSD, in 
order to recover the Einstein relation with the drift {x{t))£j^5£ — {x{t))e (here (■ ■ ■)s+5e 
denotes the average over a state where the further perturbation 5£ is applied). In what 
follows we consider that a non trivial violation of the Einstein relation happens when 
also the MSD around the drift lacks the proportionality with the drift itself. This is 
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indeed the case of superdiffusive dynamics. 

For simplicity we will refer to the case 7 = but the physical picture remains the 
same also when the case with memory is considered. In our model, by applying the 
constant field £^ > 0, we have: 

2 

{x\t)), 



2^ I v^Ti + -r, 1 

i=i ^ ^ 



4 



{Mt)]\ = {x\t))e - 



t 



(16) 
(17) 



4 (r). (r). 
In the case 2 < g < 3, namely when the the distribution Pt{t) has finite mean and 
infinite variance, the diffusion around the average position behaves asymptotically as 

'^2 ^5-5 _ t6-29 (,;2^^3-,- 



{[Sx{t)Y), 



t 



4 (r)c 

Considering for instance the case g = 5/2, by applying the matching argument to 
Eq. fITH]) . we find that the leading behaviors are 

{x\t))er^t' {Mt)]')e^t'/\ (19) 

whereas, from Eq. (fT3|) . we have that 

{x{t))eu€ - {x{t))e cc {xit))e ~ (20) 

as shown in Fig. [21 The Einstein relation is therefore violated in the out-of-equilibrium 
regime for both the MSD and MSD around the average current for all the values of the 
flight time distribution exponent (7 < 5, as summarized in Tab. [T] 





{x\t)) 






g>5 


t 


t 


t 


3<g <5 


t 


t 




2< g <3 








l<g<2 


e 




t^ 



Table 1. Columns from left to right: MSD for miperturbed dynamics; drift in the 
case of constant field £ > 0; second cumulant, i.e. {[6x{t)]'^)£ = {x^{t))£ — (a;(i)}|, in 
the presence of a constant field £. 



As already noticed in [TT] in the context of subdiffusive dynamics, such a violation 
of the Einstein relation is accompanied by an asymmetric spreading of the PDF P(x, t) 
of observing the particle in x at time t. In the case of the standard random walk, when 
a perturbation £ is applied, the P(x, t) remains Gaussian and the mean value coincides 
with the most probable one. On the contrary, in the present case, the average value of 
x{t), due to the strongly asymmetric shape, grows much faster than the most probable 
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Figure 3. Main frame: Semi-Log plot of P{x,t) at f = 10 (black histogram), and 
t = 100 (red histogram) for superdiffusive dynamics with g ~ 5/2 and constant 
external field £ = 1.5. Let us notice that for increasing time the typical (most 
frequent) value, and the mean (indicated with a vertical dashed line) of the distribution 
become increasingly different for growing times. The increasing average value of the 
distribution is due to the arising of a power-law tail for P{x,t). Inset: Semi-Log plot 
of the velocity PDF P{v). 



value. The tail of P{x, t) is reported in Fig. [3] (main frame). The stationary distribution 
of the velocities is also asymmetric, and with a power-law tail, as shown in the inset of 
Fig. [31 The study of the velocity distribution in the diffusive case can be found in [19] , 
while a careful study of the behavior of higher moments and of the scaling properties of 
P(x, t), in similar models in the absence of external field, can be found in [16] and [21] . 
respectively. 

The study of the Einstein relation on a state with non-zero current induced by a 
constant field 8 allows us to show that there is some "anomaly" in the dynamics also 
when the exponent of the power law distribution of times is g > 3. More precisely, 
when 3 < g < 5, at equilibrium, i.e. in the absence of current, a fluctuation-dissipation 
experiment would not show any anomaly in the dynamics, because {[6x{t)]'^) ~ {x(t))s. 
On the other hand, the same experiment done out of equilibrium, i.e. comparing the 
MSD around the drift with the drift itself, shows an evident violation 

mmi . (21) 

4. Conclusions 



The present study concerns the comparison of the drift of a probe particle in the presence 
of an external field with its MSD, measured either in the presence or in the absence of the 
external perturbation. The outcome of our study is that, while for a Levy walk process 
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the equilibrium fluctuation-dissipation relation is always valid [S], the situation is very 
different in the out- of- equilibrium case, where we have a breaking of the Einstein relation. 
In particular, we find that the standard fluctuation-dissipation relation can be recovered 
by replacing the MSD with the diffusion around the average current {x'^it))^ — 
as it happens for standard random walks, only for values of the exponent g > 5. This 
is a new an unexpected result, since already for g > 3 we observe a simple diffusive 
dynamics at equilibrium, even if the P{x, t) is non-Gaussian. The non-Gaussian nature 
of the distribution of displacements emerges for 3 < g < 5 only through a "response 
experiment" in presence of some currents. 

As shown in Fig. [3l the violation of the fluctuation-dissipation relation for g < 5 
comes together with a strongly asymmetric shape of the P{x, t) for large times. Let us 
note that such a mechanism, namely a transport mechanism which do not correspond 
to a uniform shifts of the P{x,t), is peculiar not only of this superdiffusive model, but 
has already been observed in the context of subdiffusive dynamics 
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